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Abstract 

We investigate the influence of the adjustment procedure and the set of measured 
observables on the properties and predictive power of relativistic self-consistent 
mean-field models for the nuclear ground state. These studies are performed with the 
point-coupling variant of the relativistic mean-field model. We recommend optimal 
adjustment algorithms for the general two-part problem and we identify various 
trends and dependencies as well as deficiencies of current models. Consequences for 
model improvements are presented. 
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1 Introduction 



There are many experimental data on the nuclear ground state. They include 
the nuclear mass, form factor, binding energy, single-particle states, half lives, 
etc. Moreover, in the future we can expect a dramatic increase in nuclear 
structure data due to new facilities that are ramping up or are envisioned 
as, for example, the future GSI facility and RIA. Our ability to understand 
the structure of the nuclear ground state is strongly connected to formulating 
nuclear models that deliver simultaneously a quantitative description of all 
measured ground-state observables. Then, present (and future) experimental 
data on nuclear ground states can be given a more realistic and complete phys- 
ical interpretation. Observables related to the mass, densities and geometries 
of the nucleus as well as specific nuclear level data can be calculated with 
microscopic mean-field models. 
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Modern self-consistent mean-field models like the Skyrme-Hartree-Fock (SHF) 
approach (for reviews see [1,2]) or the relativistic mean-field (RMF) model 
(with finite-range mesons, RMF-FR, or point-couplings, RMF-PC), see Refs. 
[2,3,4,5,6,7,8,9] have reached a mature status in nuclear structure physics. 
They are phenomenology based and are designed to allow a description of 
nuclear ground states throughout the nuclear chart, with very light systems 
(A < 16) being the exception. With few adjustable coupling constants (6- 
12) they deliver a plethora of microscopic information for a given nucleus. 
Their modern reinterpretation as Kohn-Sham schemes in the spirit of density 
functional theory [10] explains their success. They can be viewed as effective 
field theories for nucleonic degrees of freedom below an energy scale of A 
1 GeV. Short distance physics is absorbed into the various (contact-) terms 
and coupling constants. 

Independent of the interpretation of these models, each term in the Lagrangian 
or energy functional introduces a parameter or coupling constant that (yet) 
cannot be determined from more fundamental theories. Thus, the models re- 
quire parameter adjustment. A good functional form with a bad fit has small 
predictive power, thus the adjustment procedure is critical. For the description 
of exotic nuclei and reliable extrapolations to the drip-lines, a well adjusted 
mean-field is as important as a proper treatment of the pairing correlations 
and/or a possible beyond-mean-field treatment of ground-state correlations 
and excited states. 

Even though the models have reached high accuracy for various observables 
(we will discuss this in more detail later), they still exhibit significant short- 
comings and puzzles. Some of these ) the asymmetry energy of 04 ~ 
38 MeV in RMF models appears to be too large compared to other accepted 
empirical values of a 4 m 32 MeV, b) compared to SHF, the saturation density 
in RMF is smaller, the energy at saturation larger, c) there are erroneous mass 
and isovector trends in binding energies along isotopic and isotonic chains, d) 
surface thicknesses are too small (up to 5 %), e) the compressibility K seems 
to be determined in the RMF-PC variant (K « 260 MeV), while RMF-FR 
forces with values between 172 — 280 MeV exhibit roughly comparable predic- 
tive power for ground-state properties, f ) axial fission barriers for actinides and 
superheavy nuclei tend to be systematically lower for RMF compared to SHF, 
g) within RMF the neutron matter equation of state is too stiff for high den- 
sities, h) pseudo-spin splittings in RMF are too large [11], and i) the available 
types of ground-state observables are (yet) not sufficient to uniquely deter- 
mine all of the coupling constants. Note that the asymmetry energy of RMF 
models with density-dependent coupling constants is about ~ 33 MeV [12]. 
For a more detailed description of these short-comings, see Refs. [2,3,8,13]. 

In this paper we study the influence of certain methods of parameter adjust- 
ment on the resulting model, and its sensitivity to selection of observables 
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and weights. The consequences and issues related to the adjustment of phe- 
nomenological mean-field models are valid for a variety of other models and 
applications. 

We first describe the model with which we perform our studies. We then 
discuss different adjustment protocols. This is followed by specification of the 
possible sets of measured observables and their weights. The influence of these 
two topics on the resulting effective forces is then presented and discussed. 



2 Description of the Model 

For our studies we employ the point-coupling variant of the RMF model. Pre- 
vious studies [6,7,8] have shown that - with differences in detail - this model 
performs very similar to the conventional RMF-FR approach with finite-range 
fields. Yet, except for density profiles, which in the RMF-PC are smoother and 
resemble the ones obtained in SHF calculations, little sensitivity of observ- 
ables with respect to mesons versus point couplings and derivatives has been 
found. In Ref. [2], the performance of various self-consistent models, includ- 
ing the point-coupling approach, are discussed. Our current point-coupling 
Lagrangian is given by 

C = £ frcc + £ 4f + £ hot + £ dcr + £ cm , where 

£frCC _ ^([ry^Qfi _ m)lp, 

C* = -\ a s (^)(^) - | av(fy^)$7^) 

-\ a TS (V^) • ($r0) - \ a T v(^7^) • O^W), ^ 
£hot = _|/3 S (^)3 _ I 7s (^) 4 - ±7v[(V^)(^)] 2 , 

£em = - eA ^[(l - r 3 )/2] 7 ^ - \F^F^. 

We employ BCS+5-force pairing with a smooth cutoff given by a Fermi func- 
tion in the single-particle energies [14]. In each adjustment, the pairing strengths 
are adjusted simultaneously with the mean-field parameters. There are 11 cou- 
pling constants in Eq. (1) and a total of 13 coupling constants - including the 
pairing strengths for protons and neutrons - to be determined. This model has 
been employed to adjust the RMF-PC force PC-F1 [8]. Based on experience, 
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it contains a minimal number of terms necessary for a description of nuclear 
ground states comparable to established RMF-FR models. 

We note that only the direct Coulomb term is employed. We have studied 
the Coulomb exchange in Slater approximation (as it is used in SHF calcu- 
lations) and have learned that its effect constitutes a rather small correction 
for observables such as binding energies and radii. Its effect can be balanced 
by a rearrangement (readjustment) of the mean-field coupling constants. This 
does not, however, exclude the necessity for more general Coulomb exchange- 
correlation development in future models. 



3 Adjustment Protocols and Strategies 

The adjustment of mean-field models is a non-straightforward and biased en- 
deavor. First of all, there is no canonical way to do it. One has the freedom 
of choosing various kinds of observables and nuclei and of choosing the vari- 
ous uncertainties that enter the x 2 f° r the least-squares adjustment. Whereas 
some choices are obvious as, for example, binding energy as the most accu- 
rately measured ground-state observable, additional observables are not as 
obvious. Possible choices include rms radii, diffraction radii, surface thick- 
nesses, neutron radii, spin-orbit splittings, 'pseudo observables' like nuclear 
matter properties, and more. Each parameter adjustment is a balancing act 
between different types of observables, where the balance is influenced by the 
properties of the model and the weights for the various observable types. As we 
will see, current mean-field models exhibit a competition between the binding 
energy and the form factor-related observables in the fit. The nuclei chosen are 
usually spherical nuclei. Adjustments for deformed nuclei are more difficult. 

The principle problem in the adjustment of the model parameters arises from 
the complicated structure of the x 2 hyper-surface. It usually exhibits several 
separated local minima and its complicated structure results from its dimen- 
sion of D 7n 6 — 12 and the nonlinear dependence of the observables on the 
coupling constants. Often, different minima differ by the ability of the model 
to describe certain kinds of observables (energy, radii, ...) better than others. 
Also, we have sometimes found minima with quite low x 2 values and a good 
description of the fit data, but with less predictive power for the same ob- 
servables in other nuclei. We interpret these cases as an indication that our 
interaction Lagrangian £ is not (yet) complete. 

Changing the adjustment protocol, i.e., the choice of observables (see next sec- 
tion for a detailed description) changes the parameters, but this solution will 
still remain close to the point where it started, or, put differently, the relative 
minimum will still be the same when one uses down-hill methods which - in 



4 



each step - try to lower \ 2 ■ On the other hand, Monte-Carlo methods like Sim- 
ulated Annealing (SA) or Genetic Algorithms (GA) allow one to leave a certain 
region in parameter space and explore the surrounding region. Simulated An- 
nealing explores the x 2 landscape by accepting - with a certain probability 
that depends on an external temperature - parameter configurations that in- 
crease x 2 - Thus, local minima can be left and ridges can be overcome. With 
decreasing temperature, the probability for accepting higher x 2 configurations 
decreases, so that the algorithm slowly settles down in a minimum. Genetic 
algorithms obtain their ability to roam the parameter space by recombining 
different parameter vectors, their mutation and then selection of the best so- 
lutions in terms of fitness. GA can also be combined with local optimization 
with down-hill algorithms such as the ones described below. Our experience 
to date is that for new forces, a Monte-Carlo type search is mandatory for 
a first exploration of parameter space. Once several regions around minima 
have been identified, down-hill algorithms can yield the minimal x 2 values. 
The general optimization is thus a two-part problem. Since in the present 
study we are concerned with variations of adjustments starting in the same 
local minimum, we have employed only down-hill methods here. 

A variety of commercial and non-commercial optimization routines is avail- 
able on the market. Examples for downhill algorithms are simple grid searches, 
Powell's conjugate-directions method, and Levenberg - Marquardt [15,16]. The 
grid-search method, which in every iteration changes each parameter sep- 
arately before comparing the new x 2 value with the previous one is quite 
inefficient due to large parameter correlations in these models (for discus- 
sion of these correlations see Subsection 5.1). A better approach is Powell's 
method, which starts out like a grid search but adjusts the directions in pa- 
rameter space to optimize the down-hill motion. Our experience is, however, 
that in the hilly x 2 hyper-surfaces this method often gets stuck and progresses 
very slowly. Levenberg-Marquardt combines the simple gradient descent with 
a parabolic expansion of the x 2 hyper-surface around the minimum. It dy- 
namically switches between these two scenarios: far from the minimum, the 
gradient descent is preferable. Closer to the minimum, where the gradient de- 
scent would only asymptotically reach the bottom, the parabolic expansion 
is used iteratively. An additional method, denoted as trial step, uses the in- 
formation of already obtained x 2 values for a parabolic interpolation between 
them. The minimum of the parabola is the new guess for x 2 ■ 

Based on previous [8] and present experience [17], Levenberg-Marquardt plus 
the trial step in parameter space [16] (henceforth denoted as Bevington curved 
step) has been chosen here. But this algorithm needs to be restarted several 
times to reach the minimum in the x 2 hyper- surface. 

Since the x 2 landscape of these models appears to resemble a hilly landscape, 
a combination of simulated annealing and diffusion Monte-Carlo [18] has been 
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tested in [8]. This method found the region containing the parameter set PC- 
Fl, and Bevington curved step led to the actual minimum containing PC-F1. 
Algorithms and techniques such as SA, GA, Evolution Strategies, and Neural 
Networks for the optimization of x 2 might become more important with the 
increasing complexity of nuclear density functionals. They could be used for 
the optimization of the coupling constants and/or the optimization of basis 
terms in the expansion of the energy functional. 

The task of parameter adjustment is simplified if the adjustment is initialized 
with realistic parameter values, i.e., good first guesses. These initial values 
could be obtained from or related to more fundamental theories like QCD. In 
the future, as the number of terms and observables increase, such information 
will become increasingly important. One guidance has emerged already: naive 
dimensional analysis (NDA) or QCD scaling, which examines scaled coupling 
constants obtained from the original coupling constants by scaling them with 
the pion decay constant, f n = 93 MeV, and the QCD mass scale, A ps 600 — 
1000 MeV [19,20,21,22]. These scaled and dimensionless coupling constants 
are called natural if they are of order 1. So far, QCD scaling has worked well 
for recent sets of coupling constants with excellent predictive power [8,9,22,23]. 
We expect that QCD scaling will provide further useful guidance for upcoming 
more detailed model development. 

A further connection between relativistic point-coupling models and QCD sum 
rules, the chiral condensate, and pion dynamics was drawn in Refs. [24,25,26] 
where the approximate values of the coupling constants are related to the phys- 
ical properties of low-momentum QCD. Yet, due to truncated model spaces 
and other approximations, fine tuning will probably always be necessary to 
obtain optimal predictive power. 



4 Sets of Measured Observables and Weights 

The sets of observables and weights used for the adjustment of the RMF-FR 
force NL-Z2 [27] and the point-coupling force PC-F1 are shown in Table 1. 
These weights are based upon expectation and experience and constitute an 
average compromise between the actual measurement uncertainties and the 
requirement that the type of observable have a non-negligible influence in the 
X 2 minimization or adjustment procedure. The 0.2% weight on the experimen- 
tal masses is greater than the corresponding experimental uncertainties of 1.2 
x 10~ 3 % (average) and the 0.5 % weight on the experimental mean square 
charge radii is greater than the corresponding experimental uncertainties of 
0.17 % (average). 

Due to this choice of weights, the absolute value of x 2 or X 2 /dof has only a 
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relative meaning. However, the comparison of different x 2 values from various 
fits still indicates their relative goodness with respect to each other. There are 
further indications for the goodness of the fit, though. The x 2 /dof of PC-F1, 
for example, has a value of 2.75, indicating that the model does not fulfill our 
expectations of performance reflected in our choice of weights. Furthermore, 
the correlated errors of the coupling constants, which also depend on our 
choice of weights in the total x 2 i ^ill deliver relative information on how 
well the various parameters have been determined. Due to strong parameter 
correlations (see Subsection 5.1) these correlated errors range from 5-70 % 
and, in the worst cases, are larger than the value of the coupling constant, 
indicating that a readjustment without that parameter could, in principle, 
reach the same accuracy (see Ref. [8] for more details). 

We have separately reduced the weights of the observables entering our total 
X 2 , keeping the remaining weights fixed, to test the sensitivity of x 2 with 
respect to their choices. A 50 % reduction of the weight for binding energy 
(0.2 % -> 0.1 %), diffraction radius (0.5 % -> 0.25 %), surface thickness (1.5 
% — > 0.75 %), and rms radius (0.5 % — > 0.25 %) leads to an increase of the 
total x 2 by factors of 2.9, 1.4, 1.2, and 1.5, respectively. The original total 
X 2 of PC-F1 is 99. As expected, these factors are all greater than unity. The 
fact that they are (approximately) uniformly greater than unity indicates that 
(approximately) equal influence exists on the (final) interaction Lagrangian for 
each of the four observable types. This was the goal in the choice of weights, 
but is by no means unique and, above all, the choice of weights is primarily 
subject to the constraints of the actual (average) experimental errors in the 
four observable types. Repeating this exercise, but with a full parameter search 
for each of the four cases, should lead to to the same conclusion as long as one 
remains in the same minimum on the x 2 hyper-surface. 

The set PC-F1 has been used as a basis for the majority of investigations 
in this study. It has been modified for several runs to test the sensitivity 
of the resultant force with respect to various choices of observables as well 
as weights. The nuclei chosen are either spherically magic or doubly-magic 
systems. These systems are reasonably good mean-field nuclei, where ground- 
state correlations play a minor role and they require minimal computational 
time. Deformed nuclei increase the computational effort and require considera- 
tions of vibrational-rotational correlations, which are sometimes parametrized 
by phenomenological terms. However, these systems deliver additional infor- 
mation and hence will be considered in future relativistic mean-field model 
refinements, as has already been seen in non-relativistic Skyrme-Hartree-Fock 
mass formulas [28,29,30]. 

We have tested the sensitivity of the RMF-PC model with respect to the 
mass range of nuclei included in the set of observables shown in Table 1. The 
calculations displayed in Fig. 1 indicate a trade-off between light and heavy 
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PC-F1 : identical to above set 



PC-masses : binding energies only 

PC-asym : additional constraint on symmetry energy a sym = 34 MeV 

PC-surf : reduced weights on a = 0.025% 
Tabled 

Upper part: Observables and chosen weights of nuclei used for the Bevington curved 
step fitting procedure. Eb denotes the binding energy, rjr ms is the (charge) rms 
radius, i?diffr. denotes the diffraction radius, a the surface thickness, and A p and A n 
are the proton and neutron pairing gaps. A + indicates an observable contributing to 
the total x 2 - F° r the actual experimental values see Ref. [27]. Lower part: Variation 
of the sets for the four cases discussed in the present work. 

systems in adjustment procedures including nuclei ranging from oxygen to 
lead isotopes. The lighter the heaviest nuclei in the adjustment procedure are, 
the stronger the over-binding of the heavier species. A possible solution might 
be phenomenological terms accounting for correlations and Wigner energy 
and/or density dependent or gradient terms. Enhanced new-generation mean- 
field models that constitute a better approximation to the exact functional 
would presumably remedy this problem. A more pragmatic approach is to 
adjust forces for light and heavy systems separately, as done with TM2 (light 
systems) and TM1 (heavy systems) [31], but we do not pursue such a course 
here. Also, Hartree-Fock mass formulas focus only on nuclear masses. The 
approach here applies to 'complete forces' designed to describe all types of 
observables simultaneously. 



5 Results 

Certain problems in the adjustment procedure presumably arise from missing 
ingredients in the models. It is not yet clear how well correlation energies can 
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Fig. 1. Relative errors of binding energies for selected isotope (top) and isotone 
(bottom) chains for different subsets of nuclei from Table 1. 

be absorbed in the energy functional. Pairing correlations like T = pairing 
and pairing-mean-field couplings are usually missing. More generally, in both 
the isovector and isoscalar channels, terms might be missing that lead to 
deficiencies in, for example, the density dependence. The missing ingredients 
are absorbed in a somewhat uncontrolled manner in the coupling constants 
during the adjustment procedure. This, in turn, results in an entanglement 
of different parts of the effective interaction, which can then bias its physical 
interpretation. 



5.1 Parameter Correlations 



A well-known problem that occurs when adjusting mean-field models is the 
presence of strong parameter-correlations. Figure 2 shows as an example the 
total x 2 as a function of the RMF-PC coupling constants a s and ay, all other 
parameters having been fixed. These two coupling constants are responsible 
for the bulk part of the nuclear potential, and there exist correlations for the 
equivalent coupling constants in the RMF-FR and the SHF approach. The 
reason for these correlations is that not their individual values, but rather 
combinations of them have physical significance, namely, the sum as + ay 
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determines the isoscalar part of the nuclear potential, while their difference 
ots—otv determines the size of the spin-orbit potential. Thus the same potential 
depth can be achieved by various values of 0:5 and ay as long as their sum 
amounts to V xs 60 MeV. This explains the strong linear correlation. The 
limitations to this freedom arise from the density dependence of the isoscalar 
fields and the size of the spin-orbit potential which appear for larger deviations 
from the values in PC-F1. It is curious, though, that the linear correlations are 
not constant along that line in 0:5 — ay space, but that there is rather a ridge 
which rises and falls (at least) twice. This hints at further nonlinear effects 
associated with certain parameter combinations. These correlations have led 
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Fig. 2. Natural logarithm of the total % 2 for different combinations of the isoscalar 
coupling constants as and ay (here we use units of fm instead of MeV) . 

Rusnak and Furnstahl [9] to replace the original coupling constants by certain 
linear combinations of them. In the present study, however, we use the original 
parameters in order to compare to our previous adjustments. 

5.2 Weights and asymmetry energy 

Here, we study isotope and isotone chains and the impact of a variation in 
weights as well as a constraint on the value of a 4 (see Fig. 3). In these isotope 
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Fig. 3. Relative errors of binding energies for selected isotope (top) and isotone 
(bottom) chains for different choices of observables and weights, and a constraint 

On (24. 

and isotone chains, wrong trends generally indicate a problem of balance be- 
tween isoscalar and isovector fields (for isotones, also the Coulomb field). For 
example: a neutron added to a nucleus feels isoscalar attraction but isovector 
repulsion; a proton added to the nucleus feels isoscalar attraction, isovector 
attraction (since there are more neutrons than protons except for N = Z or 
proton-rich nuclei) and repulsion due to the electromagnetic field. If these dif- 
ferent fields are not balanced, an erroneous trend results. So, what we see is the 
sum of wrong trends with mass number, isospin, and electric charge. Further- 
more, the distribution of errors cannot be explained by missing ground-state 
correlations. These would produce an arch between magic corner points. Thus 
the trends in errors are genuine mismatches in the mean field. 

We can learn several things from these different adjustments. When adjusting 
only the binding energies the trends in isotopes, but even more in isotones, are 
improved. Here we have an indication that energy and form factor compete 
during the adjustment. It is interesting that in the lead isotopes beyond N = 
126 the trend seems to be reversed a bit and extrapolations might be more 
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reliable than with PC-F1. 



Including a 4 in the least-squares adjustment with a fixed value of 34 MeV 
(the usual value in RMF is around 38 MeV) worsens the isotope trends. The 
isotone trends are improved but the offset increases, and the nuclei become 
more over-bound. This may be attributed to missing repulsion coming from a 
weaker isovector channel. Other adjustments with still smaller and probably 
more realistic values of a 4 (30-32 MeV) were performed, but the resulting 
forces showed poor convergence for most non-magic nuclei. 

Putting extreme emphasis on the surface thickness a leads to very poor pre- 
dictions for the energies: for isotopes, a strong trend to under-binding results 
with increasing neutron number; for isotones, a strong trend to over-binding 
with increasing proton number results. 

The isovector repulsion drives the neutrons out. RMF models deliver a too 
small surface thickness and a too large 04. We can speculate that the too 
large 04 value in RMF models (both finite range and point-coupling variants) 
results from some missing repulsion and/or wrong density dependence in the 
isoscalar channel. A smaller value of a 4 results in less repulsion for neutrons, 
which reflects itself in the trend to over-binding in the isotope chains. In this 
way we can also understand the case with strong emphasis on a: since the outer 
part of the surface thickness is governed by the neutrons, enforcing large a 
leads to a very large value of a 4 of 40.6 MeV (see Table 2). Thus, for masses, 
neutrons suffer too much repulsion (trend to under-binding in isotopes), and 
protons suffer too much isovector attraction (see the isotones). This finding 
supports the result for the constrained 04 fit. 

Figure 4 illustrates the analysis of the form factors: adjusting only to the 
binding energies (masses), the model tends to larger radii and smaller surface 
thicknesses compared to PC-F1. A similar but not so strong effect is seen in 
the fit enforcing 04 = 34 MeV. Again, compared to PC-F1, radii are larger 
and a is a bit smaller. Enforcing an extremely high accuracy on a lifts it up to 
slightly more realistic values (there is, however, still a trend to underestima- 
tion), but at the same time leads to overestimation of the radii. This supports 
the hypothesis that the isovector repulsion is replacing repulsion that should 
come from the isoscalar channel. 

In Figure 5 we show two more quantities that yield information on the sur- 
face region of the nucleus, namely the (neutron) skin r n — r p (difference of 
neutron and proton rms radii) and the halo parameter defined as = 

r^ s — ^0.6-R^jg + 3a [32]. The most striking observation here is the effect 
of the focusing on a. While the neutron skin is slightly enhanced we see a 
strong decrease of the halo. This means that while the neutrons as a whole 
move outward, the very-low density region of the nucleus shrinks. In the cases 
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Fig. 4. Form factor observables for tin and lead isotopes for different choices of 
observables and weights, and a constraint on 04. 

of adjustments to masses only and the constrained 04, the differences with PC- 
Fl are much smaller, the skins get slightly smaller while the halo increases by 
a small amount. Again, we find an anti-correlation between halo and skin. 

Figure 6 shows the effects of the isovector terms and the emphasis on a on 
the proton and neutron densities. The main effects show up for the neutrons 
while the protons remain much less affected. The neutrons react strongly to the 
absence of isovector terms, they are pulled back into the interior of the nucleus 
and their contribution to the surface of the nucleus decreases. The emphasis 
on a in the adjustment pushes them out. Thus, any correlation between the 
neutron radius and other observables is at the same time a correlation between 
the surface thickness and these same observables. 
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Fig. 5. Halo parameters (top) and neutron skins (bottom) for tin and lead isotopes 
with the observables and weights as indicated. 



5.3 Nuclear matter 



The bulk properties of nuclear matter for the observables and weights dis- 
cussed above are summarized in Table 2. The effective mass (m*/m) is insen- 
sitive to the four adjustments, which means that the absolute magnitude of 
the scalar potential does not vary much. 

The exclusion of geometric observables from the adjustment leads to an 04 
value which is in better agreement with empirical data. This is supported by 
the adjustment with strong emphasis on the the surface - 04 is driven to 40.6 
MeV. We also observed that selections of different nuclei do not influence a 4 
drastically except that for relatively light systems a 4 decreases to about 36 
MeV. 

While the typical value for the compressibility in the RMF-PC approach is 
K 7n 260 — 280 MeV, this value is lowered to probably even more realistic 
values in two adjustments, namely, the a fit and the fit only to masses with a 
percentage error. 
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Fig. 6. Proton and neutron point densities (left) and their differences (corresponding 
to ptv) for different choices of observables and weights with respect to the isovector 
channel and surface thickness a. 



Type of adjustment 


Po [fm -3 ] 


E/A [MeV] 


K [MeV] 


m* jm 


a 4 [MeV] 


PC-F1 


0.151 


16.17 


270 


0.61 


37.8 


PC-masses 


0.150 


16.03 


237 


0.61 


33.9 


PC-asym 


0.151 


16.12 


268 


0.61 


34.0 


PC-surf 


0.147 


16.28 


238 


0.62 


40.6 



Table 2 

Bulk properties of nuclear matter for the various adjustment protocols discussed in 
the text. 

6 Conclusions 

In this paper we have dealt with a less frequently discussed but neverthe- 
less quite important topic in phenomenological nuclear structure, namely the 
adjustment of the model parameters (coupling constants). We have studied 
the predictive power of the relativistic mean-field model with point-couplings 
determined from various nuclear ground-state observables. The dependencies 
on the weights of different observables in the x 2 minimization and on mass 
ranges of nuclei have been evaluated. Alternative adjustment protocols have 
been tested and issues as well as features of current adjustment algorithms 
have been discussed. 



15 



We have tested various optimization algorithms and found the Levenberg- 
Marquardt method, supplemented with a trial step in parameter space (Bev- 
ington curved step), to be the most efficient method. This is true for adjust- 
ments aiming at the nearest minimum. The x 2 hyper-surfaces of nuclear mean- 
field models, however, contain several (local) minima. Thus, for adjustments 
of new forces, Monte-Carlo methods like simulated annealing are necessary 
to explore the (unknown) parameter space. Once the various minima have 
been identified, a down-hill method can drive the parameter vector into the 
minimum. 

The main outcome of these studies is that the present RMF forces are not 
able to describe simultaneously the mass and the geometry of the nucleus 
with the desired accuracy. The problem appears to be the surface thickness 
of the nucleus. It is underestimated in all forces under investigation, and this 
underestimation is linked to the corresponding asymmetry energies. In the 
present sets, a larger and more realistic surface thickness leads to an even larger 
value of 04, which is already too high in RMF models (04 38 MeV) when 
compared with values from liquid-drop mass fits yielding 04 ~ 32 MeV. This 
finding supports earlier results on a correlation between CI4 and the neutron 
skin [33]. 

Relativistic mean-field models are aimed at describing nuclei from A > 16 
up to superheavy nuclei with A ps 300. Thus, the current mass range in our 
adjustments is from oxygen to lead isotopes. Our results show that the present 
forces constitute a compromise of describing both light to medium and heavy 
systems. Heavy systems get over-bound if the adjustment is confined to lighter 
nuclei. Correlations in lighter systems, which constitute a larger fraction of the 
binding energy compared to heavy systems, demand stronger binding leading 
to too strong binding in heavy nuclei. Hence, current adjustments have to 
live with this trade off. Possible solutions are more terms that can account 
for these missing correlations and/or more phenomenological add-ons like a 
Wigner-energy term to account for missing Z = N pairing. 

More generally, our investigations indicate that the isoscalar and isovector 
channels are somewhat entangled in current forces and that hence the isovector 
channel might account for mismatches in the isoscalar part of the interaction. 
To avoid this entanglement, more observables of isovector type are needed, 
especially from nuclei with very large isospin. 

It is clear that additional terms are needed to account for the missing physics 
in the present RMF approaches. They might - to a large extent - account for 
the missing correlations and density dependence. Alternatively, more detailed 
range might be required. A systematic exploration of various extensions in the 
framework of RMF-PC is underway [34]. 
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